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Abstract. The general conclusion of Seevinck and Larsson is that our model exploits the so-called 
coincidence-time loophole and produces sinusoidal (quantum-like) correlations but does not model the 
singlet state because it does not violate the relevant Bell inequality derived by Larsson and Gill, since in 
order to obtain the sinusoidal correlations the probability of coincidences in our model goes to zero. In this 
reply, we refute their arguments that lead to this conclusion and demonstrate that our model can reproduce 
results of photon and ion-trap experiments with frequencies of coincidences that are not in conflict with 
the observations. 
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In order to come to their conclusions, Seevinck and 
Larsson made the following statements pQ: 

— De Raedt et al. claim that their model violates the 
CHSH inequality, a claim that cannot be found in 
Ref. [2]; 

— The CHSH inequality is inappropriate for models that 
exploit the so-called coincidence-time loophole [5] and 
the appropriately modified inequality [3] is not violated 
by the model of De Raedt et al.; 

— The model of De Raedt et al. cannot reproduce all 
experimental realizations of the EPRB experiment; 

— De Raedt et al. claim that their model can repro- 
duce the coincidences of recent experimental results, 
another claim that cannot be found in Ref. [2]; 

and put our model in the context of hidden variable mod- 
els to obtain an expression for the probability of coinci- 
dences. 

In this reply, we point out once more that in our work [2] 
we did not rely on Bell's or CHSH's inequality nor on any 
generalization thereof to come to our conclusion that it 
is possible to construct an event-based simulation model 
that satisfies Einstein's criteria of local causality and real- 
ism and can reproduce the expectation values of a system 
of two 5=1/2 particles in the singlet state. In our work [2] 
we did not make any claims about these inequalities, nei- 
ther did we make any statement about coincidences in real 
experiments. We furthermore demonstrate that the calcu- 



lation by Seevinck and Larsson [T| of the probability of co- 
incidences for our model is simply wrong and we present 
results for the frequency of coincidences which compare 
rather well to the values observed in recent experiments. 

Before replying in detail to the comments of Seevinck 
and Larsson [1] on our model [2], we first want to sin- 
cerely apologize that we did not make a reference to the 
model presented in Ref. [3] which, like our model, uses 
coincidence in time as a criterion to decide which pairs of 
detection events are to be considered as stemming from 
a single two-particle system. Furthermore, to our knowl- 
edge, it is the first work to point out how the original Bell 
inequality changes when using this post-selection proce- 
dure |3J. Although in our work [2J, we did not rely on Bell's 
original inequality or on any of its generalizations [4] to 
come to our conclusions, we should have made reference to 
Ref. [3] only because of the fact that the model presented 
in Ref. [3j uses the same pair selection criterion as we use 
in our model [2]. 

As stated in paragraph 2 of our paper [5] , we consider 
the original EPRB problem, that is the construction of a 
model that satisfies Einstein's criterion of local causality 
for each pair of events and reproduces the expectation 
values of a system of two 5 = 1/2 particles in the singlet 
state. Whether or not this model leads to a violation of 
some inequality is of secondary interest. 

We do not share the point of view of Seevinck and 
Larsson p] that a system is in the "singlet state" if and 
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only if some correlation violates a certain bound and if 
the probability of coincidences does not go to zero. This 
viewpoint does not make much sense if we apply it to the 
ground state of a hydrogen molecule (a spin singlet), for 
instance. The singlet state is a concept of quantum theory. 
Unlike Seevinck and Larsson suggest in their conclusion, 
the singlet is not defined with any reference to coincidence 
counts or (generalized) Bell inequalities. 

According to quantum theory, the singlet state is com- 
pletely characterized by the single-particle and two-particle 
expectation values ia(ai) = -^2(^2) = and E(&i,a.2) = 
— ai-a2, respectively. Because quantum theory has nothing 
to say about single events [5] , it does not give us a recipe 
to compute -Ei(ai), £^(^2) an d B(ai,a2) from the record 
of single events in laboratory experiments or theoretical 
models. 

In the case of EPRB laboratory experiments with pho- 
tons [B] , coincidence in time seems to be a good criterion 
to identify particle pairs based on the time-tag data of 
single particle events, since by using this criterion results 
comparable to those expected from quantum theory can 
be obtained. Therefore we use the same criterion in our 
model [2]. However, other criteria to decide which single 
particles belong to a single two-particle system are not ex- 
cluded. The criterion depends on the experimental setup 
but quantum theory does not give any guidance to define 
a criterion. Having made a definite choice for this crite- 
rion, we can compute the single-particle and two-particle 
expectation values from the record of single events in labo- 
ratory or computer experiments and compare the outcome 
with the results from quantum theory. If and only if we 
find i?i(ai) = E 2 {a2) = and £ , (ai,a 2 ) = ai • a 2 , we 
may say that we found expectation values that correspond 
to those of a singlet state. No other criteria, like violating 
an inequality or computing a probability of coincidences 
for example, are required to come to this conclusion. Note 
that we cannot say anything more than that the expec- 
tation values correspond to those of a singlet state. For 
example, we cannot make statements such as the source 
produces singlets, since the results for the expectation val- 
ues do not only depend on the characteristics of the source 
and the detection elements but also on the measurement 
(post-processing) process. 

It is self-evident that our model is too simple to de- 
scribe, in every detail, all conceivable experimental real- 
izations of the EPRB thought experiment but it is the first 
model that satisfies Einstein's conditions of local causal- 
ity and realism and that exactly reproduces the single- 
particle and two-particle expectation values of the singlet 
state. In this respect, it may be viewed as the first real- 
ization of the EPRB thought experiment (as defined by 
EPRB), since none of the laboratory experiments of the 
EPRB experiment have shown results for the single- and 
two-particle expectation values that compare so well with 
those of quantum theory. In those experiments, conclu- 
sions are usually drawn based on the value for S max only. 
Moreover, to draw conclusions about local realist mod- 
elling of expectation values that agree with those of a sin- 
glet state, finding one such model is sufficient. Whether 



this model then fails to describe all possible laboratory 
realizations of the EPRB thought experiment becomes 
irrelevant and it remains to be seen if these laboratory 
experiments produce data that completely characterize a 
singlet state, a requirement of the EPRB thought experi- 
ment. 

Seevinck and Larsson state "we will put the model used 
by De Raedt et al. in its proper context" pQ. However, they 
failed to do so in any respect. In spite of the fact that in 
our paper, we repeatedly stress that in formulating our 
model we do not rely on concepts of probability theory, 
they seem to ignore our statements. This is most evident 
by their statement that "The local hidden variable ... is 
denoted by 5 n ,i". Since, according to Larsson a hidden- 
variable model is really a probabilistic model [7] and since 
our model is purely ontological, the concept of a hidden 
variable cannot be applied to our model as such. 

In contrast to the (repeated) statement made in Ref. pQ , 
we did not claim that the CHSH inequality (see Eq. (2) in 
Ref. pQ) is valid for our model. There is no such statement 
in our paper. In our paper, we studied the values of S max 
as a function of the time window W relative to the time- 
tag resolution r and this for several values of the model 
parameter d. We compared the results with S max — 2\/2, 
the quantum theoretical result for the singlet state and 
also the maximum value for S max that can be obtained 
for any choice of the quantum state. Although we find 
that for some model parameters 2 < S max < 4 we did 
not claim that our model violates the CHSH inequality, 
as stated by Seevinck and Larsson pQ. In fact, using ele- 
mentary algebra it follows immediately from Eqs. (3) and 
(5) of Ref. that |_E(ai,a 2 )| < 1 and that 

\E(&, c) - E{&, d) + E{b, c) + E(h, d)j < 4, (1) 

for the data generated by our computer model. Without 
any further constraints on the algorithm that generates 
the data T 2 } (see Eq. (1) in Ref. [2]), the upperbound 
(4) in Eq. (pQ) cannot be improved. In our paper [2], we use 
expression Eq. (5) (see Ref. [2]) to discuss the nature of 
the quantum state, but attach no meaning to the violation 
of some bound by our simulation data. 

Not surprisingly, also the statement "They further- 
more claim that the maximal quantum violation is pQ 
is wrong. We did not make any reference to the CHSH in- 
equality in our paper. Looking at Fig. 3 of our paper [2], 
Seevinck and Larsson should have noted that for d > 
3, our model can produce correlations that are (much) 
stronger than those of quantum theory (which in view of 
Eq. (pQ) is not a surprise). In fact, the correct statement 
(see Ref. [2]) is that our model can exhibit correlations 
that are stronger than those of quantum theory of two 
5=1/2 particles. 

As we mentioned before, we agree with Seevinck and 
Larsson that when time-coincidence is used to decide which 
pairs of detection events are to be considered as stemming 
from a single two-particle system and if one would like 
to consider a generalized Bell inequality, the relevant in- 
equality to consider would be Eq.(4) in Ref. pQ and not 
the CHSH inequality. In this modified inequality 7 is the 
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infimum of the probability of coincidence |3| . Seevinck and 
Larsson compute 7 for our model. They conclude that our 
model does not violate Eq.(4) in Ref. pQ. Moreover, al- 
though in the paper on the photon experiment [6] or in 
the paper on the ion-trap experiment [5], there is no in- 
formation about the minimum frequency of coincidences 
(the minimum is required for the application of Eq.(4) 
in Ref. pQ), Seevinck and Larsson refer to these papers 
when they cite the values of 7 = 0.05 and 7 = 1 and then 
state that our model cannot reproduce the frequencies of 
coincidences that agree with those found in these two ex- 
periments. First of all, the statement that "De Raedt et 
al. claim that their model can reproduce the coincidences 
of recent experimental results" [1] is simply wrong: There 
is no such claim in our paper [2]. Second, it is logically 
inconsistent to draw conclusions based on the compari- 
son with the ion-trap experiment for which 7 = 1 |TJ and 
third, the calculation of 7 [lj for our model is wrong, two 
statements which we will prove in what follows. 

Let us formalize the statements in Ref. pQ as proposi- 
tions (denoted by A,B,...): 

A. The probabilistic, hidden variable model using the time 
window to define coincidences [3] yields the inequality 
given in Eq.(4) of Ref. PQ (see also Ref. [5]). The up- 
perbound of this inequality is given by 6/7 — 4, where 7 
denotes the infimum of the probability of coincidences 
over all possible settings ai, &2 of the detectors. 

B. In the EPRB experiment with ions [8], 7 = 1 pQ. 

C. The probabilistic, hidden variable model of Ref. [3] ap- 
plies to the ion-trap experiment [8J and hence the ex- 
perimental data should satisfy the inequality given in 
Eq.(4) of Ref. pQ. Note that the second part of this 
statement implicitly follows from proposition B. 

D. The ion-trap experiment yields S max ~ 2.25 [SJ . Strictly 
speaking, this statement is not made in Ref. pQ, but it 
is an experimental fact and as such cannot be denied. 

Let us now apply the rules of elementary logic. 

If 7 = 1 pQ , the ion-trap experiment not only violates 
the original Bell inequality but as 6/7 — 4 = 2, it also 
violates the inequality given in Eq.(4) of Ref. pQ. Thus, 
we have 

AABACAD^C, (2) 

where A, =>■ and denote the logical "and" operation, log- 
ical implication, and logical negation, respectively. Clearly, 
Eq. ([2]) expresses a logic contradiction. If we assume that 
propositions A and D are true (as we do), then we must 
conclude that B or C or both B and C are false. In any 
case, the argument used by Seevinck and Larsson leads to 
a logical contradiction, independent of what we wrote in 
our paper [2J. 

We should not exclude the possibility (that is, we might 
accept proposition C) that the model of Ref. [3] or ours [2] 
is too simple to describe the ion-trap experiment 8J . This 
experiment uses a detection pulse during which the bright 
state of an ion scatters many photons (64 on average) [5]. 
This process may not be sampling "single-events" but is 
more likely to probe the ensemble average that is given by 



quantum theory (although the number of samples, ss 64, 
is not large). 

By trying to put our work in the context of "hidden 
variable theories", Seevinck and Larsson also made mis- 
takes in elementary algebra. Seevinck and Larsson assume 
that the probability of coincidences is given by the denom- 
inator of Eq. (6) in Ref. [2] (see Appendix A of Ref. pQ), 
from which they derive an expression for the probability of 
coincidences 7 (see Eq. (8) in Ref. pQ). However, Seevinck 
and Larsson apparently overlooked the fact that in going 
from Eq. (3) to Eq. (6) (see Ref. [2J), we take the limit 
W/Tq — t/T — > and let the number of events N in 
both the numerator and denominator go to infinity. Al- 
though the ratio remains finite, which is obvious in the 
case ai = a.2 (x±X2 = —1) where it is equal to minus one, 
the limit of the denominator may not exist and in fact, it 
diverges if ai = a.2- This is not a problem of our model: 
This divergence merely signals that one has to be care- 
ful in taking the limits. The mathematical derivation in 
Appendix A of Ref. pQ is simply incorrect. 

Nevertheless, Seevinck and Larsson raise an interest- 
ing question about the role of the frequency (not proba- 
bility) of coincidences in our model. For nonzero time-tag 
resolution r and time window W > r, the frequency of 
coincidences in our simulation model is given by 

1 N 

r = 1 j J E&(w-\t n , 1 -t n , 2 \), (3) 

n=l 

a well-defined quantity in our simulation model that is 
easy to compute numerically. Notice that r is a func- 
tion of ai and &2 and that < -T < 1. Assuming that 
the results we obtain by using pseudo-random numbers 
can be described by a probabilistic model, we expect that 
7 = min ai !32 r with probability one if N is sufficiently 
large. With these additional assumptions, not only the in- 
equality Eq. (pQ) holds but also the inequality given by 
Eq. (4) of Ref. pQ holds. 

In our model, there are four free parameters, namely 
the time window W/t, the maximum time delay Tq/t, 
the time-tag exponent d and the number of events N. 
For d = 3, N -> 00 and in the limit W/T = t/T -> 
0, our model reproduces exactly, the expression for the 
two-particle expectation value of a quantum system in the 
singlet state [5J. 

For d = 3, W = T, T /t = 1000 and N = 10 6 (the 
results reported in this paper do not change if N > 5 x 
10 5 ), we find that min aia2 r s» 1.27W/Tq, in concert with 
the rigorous result (for d = 3 and W = r) 

min <^ lim r } = -— » 1.27—. (4) 

ai,aa [w/T o ^0 J 7T T T 

Thus, in the regime W/T — t/T — > 0, we find that the 
minimum frequency of coincidences is proportional to the 
width of the time bins, as it should be. 

Next, we consider the possibility of fitting the results 
of our model to the experimental data of an EPRB exper- 
iment with photons ,6 a and ions [8] . According to Seevinck 
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Fig. 1. (color online) The frequency of coincidences f as a 
function of 9 = arccos(ai ■ a.2) for parameters r/Tb, W/To and 
d chosen such (see text) that the simulation model reproduces 
the result for a singlet state, S m „ = 2.83 (solid line, red), and 
the values of Smax = 2.25 (dashed line, black) and Smax = 2.73 
(dotted line, blue), as obtained from experiments with ions [8] 
and with photons [6], respectively. 
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Fig. 2. (color online) Simulation results of the two-particle 
correlation _E(ai,a2) as a function of = arccos(ai ■ a2) for 
the model parameters that yield Smax = 2.25 (squares, black), 
Smax = 2.73 (stars, blue), and Smax = 2.83 (bullets, red), 
respectively. The solid line (black) is the result (_E(ai,a2) = 
ai ■ a2) for the singlet state. 



and Larsson pQ , our model cannot reproduce these exper- 
imental data. For simplicity, we set Tq/t — 1000, d = 3 
and take N — 10 6 . Then, there is one free parameter left, 
namely the (dimensionless) time- window W/t > 1. The 
fitting procedure consists of changing WJ r such that the 
value of S m ax = max ai b )C ,d «5(a, b, c, d) agrees with values 
cited in Refs. [6"[[8]. 

In Fig. [TJ we present our results for the frequency of 
coincides for the values of W/t — 285, W/t — 16, and 
W/t — 1, for which our model yields S max — 2.25 [8], 
Smax = 2.73 [6], and S moi = 2.83 (singlet state), respec- 
tively. 

From Fig. [TJ we see that in order to reproduce the ion- 
trap result [8], the frequency of coincidences r > 0.52 is 
quite large. It is important to recognize that with four free 
parameters at our disposal, it is easy to reproduce almost 
any number for S m ax, as long as it is between zero and 
four. For instance, we find the same value of S m ax — 2.25 
for W — t and T /t = 1.025 but then T > 0.87. In any 
case, these results refute the statement in the Comment 
that our model cannot reproduce the experimental result 
Smax — 2.25 of the ion-trap experiment [5] with a nonzero 
value of r. Fitting our model (for d = 3 and T q /t = 1000) 
to Smax = 2.73 and S ma x — 2.83 yields F > 0.0377 and 
r > 0.00127, respectively. 

An analysis of experimental data for an EPRB ex- 
periment with photons [9] yields r w 0.01 (the value of 
7 « 0.05 cited in Ref. [TJ is the total frequency of coinci- 
dences, that is the sum over four experiments, and not the 
infimum over all possible experiments, as required for the 
application of the inequality given in Eq. (4) of Ref. [TJ). 
Thus, for the same value of S m ax , our model yields a value 
of r that is larger (r — 0.0377) than the value that can 
be extracted from experimental data for an EPRB exper- 
iment with photons [9]. 



As our model is flexible enough to yield for S m ax any 
number between zero and four with reasonable values of 
the model parameters, it is of interest to study how the 
correlation _E(ai,a2) deviates from the result £7(ai,a2) = 
ai • a2 of a system in the singlet state as we fit the values 
of S. m ax to the experimental results. 

In Fig. [2j we show the simulation results for the same 
three cases S ma x = 2.25, 2.73, 2.83. From Fig. [2j it is clear 
that the simulation data that yields S max — 2.25, 2.73 
cannot be described by a single sinusoidal function, but 
for Smax > 2.73 the deviations from a single sinusoidal 
are small and it remains to be seen if experiments can 
resolve such small differences. 

As is evident from Fig. 3 in Ref. [2], for d > 3 our 
model yields the value for the singlet state S ma x = 2\f2 
without having to consider the limit W/T = t/T — > 0. 
Thus, in order for an experiment and a model of the type 
considered in our paper to reproduce all the features of 
a quantum system of two S — 1/2 particles in the sin- 
glet state, it is not sufficient to show that it can yield 
Smax = 2\f2 for some choice of the parameters. As men- 
tioned before, the singlet state is completely characterized 
by the single and two-particle expectation values. Hence, 
in order to make a comparison with the singlet state, it is 
necessary to measure or compute these two quantities. 

Finally, the statement in our paper [5] that "Our work 
suggests that it is possible to construct event-based sim- 
ulation models that satisfies Einstein's criteria of local 
causality and realism and can reproduce the expectation 
values calculated by quantum theory |10illlill2ill3l[T4] " 
should not be taken out of the context as Seevinck and 
Larsson did by omitting the references. In fact, what we 
have shown in the work that we refer to is that it is possi- 
ble to perform an event-based simulation, satisfying Ein- 
stein's criteria of local causality, of a universal quantum 
computer [13) , which according to the theory of quantum 
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computation should suffice to simulate any quantum sys- 
tem [T5] . 

In conclusion, the purpose of our work is to construct 
an event-based simulation model, satisfying Einstein's cri- 
teria of local causality and realism, that produces the 
quantum correlations of the singlet state [2] . As we have 
shown in Ref. [2] we succeeded, to our knowledge for the 
first time, in constructing such a model. Our conclusion 
that we find results that are indistinguishable from those 
of a singlet state is based on the fact that the calculated 
single-particle averages and two-particle correlation func- 
tion agree with the well-known results £i(ai) = £^2(^2) = 
and -E , (a 1 ,a2) = — ai ■ a2 for a system of two S — 1/2 
particles in the singlet state. 

We have also demonstrated in Ref. [2] and in this re- 
ply that knowing S max , a quantity derived from the two- 
particle correlation function, does not suffice to draw any 
conclusion about the observation of a singlet (-like) state. 
We also demonstrate in this reply, that our model can not 
only produce the results from quantum theory for a sys- 
tem of two S — 1/2 particles in the singlet state but that it 
can also be applied to EPRB laboratory experiments with 
photons and ions and give results for S max and the fre- 
quency of coincidences that are comparable to the values 
extracted from these experiments. 
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